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Abstract 
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1 Introduction and the main results 



The classical invariant ensembles of random matrices are given by probability measures of the 

form e -NfSTrV{H)/2 where N jg the 

size of the matrix H and V is a real valued potential. The 
parameter j3 — 1, 2, 4 is determined by the symmetry type of the matrix, corresponding respectively 
to the classical orthogonal, unitary or symplectic ensemble. Let A = (Ai, A2, . . . , Ajv) 6 Sjv be 
the eigenvalues of H in increasing order, where Sjj C M. N denotes the simplex determined by 
Ai < A2 < . . • < Xn- It is well- known that the probability distribution of the ordered eigenvalues 
on Ejy is given by 

M W =/i W(dA) = ^ e -W)dA, H(X) = f; \v{X k ) - I Y. MAi-A,). (1.1) 

N k=l l<i<j<N 

For non-classical values of j3 > 0, i.e., j3 $ {1,2,4}, one can still consider the measure (1.1) on 
Ey, but in general there is no simple natural matrix model producing this measure except for 
the Gaussian case, V(x) = x 2 , which corresponds to a tri-diagonal random matrix [13,26]. We 
will view n = ^ N > as a Gibbs measure of particles in R with a logarithmic interaction, where the 
parameter (3 > is interpreted as the inverse temperature. We will refer to the variables Xj as 
particles or points and the system is called log-gas or general /3-ensemblc. 

The universality conjecture asserts that the eigenvalue gap distributions in the bulk depend only 
on f3 and are independent of the potential V. For classical ensembles, the eigenvalue correlation 
functions can be explicitly expressed in terms of polynomials orthogonal to the measure e~@ v ( x ^ 2 . 
Thus the analysis of the correlation functions relies heavily on the asymptotic properties of the 
corresponding orthogonal polynomials. This approach, initiated by Dyson, Gaudin and Mehta (see 
[21] for a review) was the starting point for all results on classical universality Precise analysis on 
orthogonal polynomials for general class of weight functions was made possible by the Riemann- 
Hilbcrt approach [5,11,12]. There are also methods independent of the Riemann-Hilbert approach, 
see, e.g., [20,22,23]. The universality for j3 — 2 was proved for very general potential. For ,9 = 1,4 
[10, 18, 25] it was proved for analytic V with some additional conditions. A summary of recent 
developments can be found in [2,9,10,24]. 

For non-classical values of /3, i.e., /3 $ {1, 2, 4}, there is no simple expression of the correlation 
functions in terms of orthogonal polynomials. In [6], we initiated a new approach to prove bulk 
universality for all ft > and strictly convex V. The method was based on estimating correlation 
functions by local Dirichlet form and the main ingredients consist of the following two steps: 

Step 1. Rigidity of eigenvalues. This establishes that the location of the eigenvalues are not too 
far from their classical locations determined by the equilibrium density p(s). 

Step 2. Uniqueness of local Gibbs measures with logarithmic interactions. With the precision of 
eigenvalue location estimates from Step 1 as an input, the eigenvalue gap distributions are shown to 
be given by the corresponding Gaussian ones. (We will take the uniqueness of the gap distributions 
as our definition of the uniqueness of Gibbs state.) 

Our goal is to extend this result to the non-convex case. It was emphasized in [6] that the 
convexity of the potential V was used only in Step 1. So in order to apply this method, it suffices 
to prove the rigidity estimate which we now introduce. 

We will assume that the potential V is real analytic function in R such that its second derivative 
is bounded below, i.e. we have 

inf V"(x) > -2W (1.2) 

16E 
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for some constant W > 0, and 

V(x)>(2 + a)ln(l + |x|), (1.3) 

for some a > 0, if \x\ is large enough. It is known [7] that under these (in fact, even weaker) 
conditions the measure is normalizable, < oo. Moreover, the averaged density of the empirical 
spectral measure, defined as 

p[ N \x) = p{ N ^(X) :=E, w lf>(A-A,) 

converges weakly to a continuous function p, the equilibrium density, with compact support. We 
additionally assume that p(s) is supported on a single interval [A, B], and that V is regular in the 
sense of [19]. We recall that V is regular if its equilibrium density p is positive on (A,B) and 
vanishes like a square root 1 at each of the endpoints of [A, B], that is 

p{t) = s A Vt^A[l + 0(t— A)) , t^r A + , (1.4) 



p{t) = s B VB-t (1 + (B -t)), t->B 

for some constants sa, sb > 0. 

In this pap 
and defined by 



In this paper, we are interested in the usual n-point correlation functions, generalizing p^\ 



pW(x 1 ,...,x n )= p(x)dx n+1 . . .dx N , (1.5) 



where p is the symmetrized version of p given in (1.1) but defined on K instead of the simplex 
Eat: 

where A (<t) = (A (t( i ) , . . . , X a (N)), with A ct(1 ) < • ■ ■ < X a {N)- 

In the following, we omit the superscript N and we will write p for p( N >. We will use P M and 
E M to denote the probability and the expectation with respect to p. Let the classical position 
be defined by 

k 

p(s)ds = — . (1.6) 



N 

Finally, we introduce the notation \p, q\ = [p, gjflZ for any real numbers p < q. 

It is known that the particles are rigid, i.e. they cannot be far from their classical locations 2 : 
for any e > there are positive constants c\, C2 such that, for all N > 1, 

P„(3fc€ [1,JV] I |A fc -7fc| >e) < Cl e- c * N . (1.7) 



1 This is not a strong constraint: [19] proves that the regular potentials V are a dense and open subset of the 
potentials for the topology induced by the distance 

l + ll^ (j, -^0)|| L » [ _ feifc] j^J l + |G fc (V)-Gfc(W0r 

where Gfc(V) = inf^^j. V(a;)/log 

2 For eigenvalues in the bulk, (1.7) follows from the large deviations for the empirical spectral measure with speed 
N 2 [2,4], and for the extreme eigenvalues the large deviations principle with speed N was proved in [2], Theorem 
2.6.6, up to a condition on the partition function that follows from Theorem 1 (iii) in [25]. 
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The main technical result of this paper is to prove that rigidity holds for the measure p at the 
optimal scale 1/N in the bulk in the following sense. This theorem extends our rigidity result in 
[6] to non-convex potential V. 

Theorem 1.1 (Rigidity estimate in the bulk) Let V be real analytic, regular with equilibrium 
density supported on a single interval [A,B], and satisfy (1-2), (1.3). Take any a > and e > 0. 
Then there are constants S,c\,C2 > such that for any N > 1 and k £ [aiV, (1 — a)N\, 

P M (|A* - 7fc| > ^~ 1+£ ) < c ie - C2N \ (1.8) 

Our main result on the universality is the following theorem: 

Theorem 1.2 (Bulk universality) Let V be real analytic, regular with equilibrium density sup- 
ported on a single interval [A,B], and satisfy (1-2), (1-3). Then for any j3 > the bulk universality 
holds for the [3-ensemble p = pp,v- More precisely, for any E € (A, B) and \E'\ < 2, for any 
smooth test functions O with compact support and for any < k < \, we have, with s := N~ 1+k , 
that 



lim / d«i • • • da n 0{a\, . . . , a T{ 

JV-s-oo J 



E+s 

E-s 



dx 1 (Ar) / oi a„ \ 

2sp(E) nPn \ X Np(E) 7 '" ,X Np(E)J 



E +s dx 1 (AT) / ol\ a r , 



(N) ( ai , 



,_ s 2sp sc {E'Y PGauss -"V N PBe (E r y'"' N Psc (E>) 



0. 



Here p sc (E) = J--\f£ — E 2 is the Wigner semicircle law and p[^2 uss n are the correlation functions 
of the Gaussian (3-ensemble, i.e. with V(x) = x 2 . 

Theorem 1.2 follows immediately from the rigidity estimates, (1.7), (1.8), and the uniqueness 
of local Gibbs measure, i.e., Theorem 2.1 and Corollary 2.2 in [6]. We note that the proof of the 
latter results in Section 4 of [6] uses only the rigidity estimate, given in Theorem 3.1 of [6], as an 
input. Once the rigidity estimate is proven, the rest of the argument is identical and we will not 
repeat it here. 

The rest of this paper is devoted to the proof of Theorem 1.1. After some initial estimates 
concerning the large deviations regime and global smooth linear statistics (Section 2), the proof 
consists in the following steps. First we compare p to some convexified measures v (Section 3); 
the Hamiltonian of v differs from that of p mainly by some properly chosen linear statistics 
of the A-s's, allowing % v to be convex. Despite this change in convexity, we will prove that the 
two measures p and v have the same subexponentially small probability events. This step is the 
main extra ingredient allowing one to generalize the rigidity estimate obtained in [6]. Then by a 
self-improving method, this measure v (together with p) is proved to have rigidity till the optimal 
scale, thanks to comparisons with locally constrained versions of v (Section 4). 



2 Preliminary results 

2.1 Equilibrium measure, large deviations 

For analytic potential V satisfying the asymptotic growth condition (1.3), the equilibrium measure 
p{s)ds associated with (p^ n ^)n>o can be defined as the unique minimizer (in the set of probability 
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measures on E endowed with the weak topology) of the functional 

I{v)= ( V{t)dv{t)- [[ log |t — s\dv(s)dv(t) 



if J V(t)dv(i) < oo, and I(y) = oo otherwise. Moreover, if one assumes that p is supported on a 
single interval [A, B] and regular in the sense of the previous section, p has the following properties: 



(a) This equilibrium measure satisfies 



iV( t , - / p.!) 



for any i <G (A, £?). 
(b) For any f € [A-B], 

p(t) dt = ?-r(t)y/(t-A)(B-t)l [AtB] dt, (2.2) 
where r can be extended into an analytic function in C satisfying 

r(z) = -1 f K^Tft) d * (2 . 3 ) 

2tt J A z-t y/(t-A)(B-t) 

In order to have the density supported strictly in a compact interval, for given re > 0, define 
the following variant of p^ N ^ conditioned to have all particles in [A — re, B + «]: 



iV 

/i W«)( dA)= \Xi-Xj\ P l[e- N % nXk h Xke[A ^ B+K] d\ 1 ...d\ N . (2.4) 



l<i<j<N fe=l 

In this paper we will choose « to be small. This choice differs from [6] where, instead of 
[A — K, B + k], we restricted the particles to [-R, R] for a very large R. The smaller interval is 
needed here because we need r to be positive on the support of ^( Ar ^ K ) in the proof of Lemma 2.2. 
Unlike in the case of convex V where r is known to have no real zero at all, for the non-convex 
regular case we only know that r is nonzero in the interval [A, B\. By continuity, it is also nonzero 
in [A — k, B + re] for some small k. 

Let p\^'^ denote the correlation functions of the measure p( N ' K \ Then Lemma 1 in [7] states 
that under condition (1.3), for some large enough re there exists some c > 0, depending only on V, 
such that for any xi, . . . , Xk £ [A — re, B + re], we have 

Pk (xi,---,x k )-p[ N) (x 1 ,...,x k ) <p k N ^\xi,...,x k )e- cN , (2.5) 

and for x\, . . . , xj £ [A — re, B + re], Xj+\, . . . , x k S [A — re, B + re], 

pf ) (x 1 ,...,x fc )< e -^^ 1 1 °sl-l. (2.6) 

The estimates (2.5) and (2.6) actually also hold for arbitrarily small fixed re > thanks to the 
large deviations estimates (1.7). 
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2.2 Linear statistics 



The following lemma was essentially proven in [25] (for the variance of linear statistics). 



Lemma 2.1 For any {unction (f> with \\<j. 
depending only on V and (f> ( one can choose c 
N > 1 and s > 0, 



+ H^'lloo + ||0"||oo < oo, there is a constant c > 



O(||^||oo + 1 1 <A' | |oo + ||<^"||oo)j such that, for any 



N 

S24>(Xi)-N / p(u)<j>(u)di 



> s < e 



-cs/ log N 



Proof. Without loss of generality, we can assume that <f> is compactly supported (thanks to 
large deviation estimates such as (2.6)). We know from Shcherbina, equation (2.22) in [25], that 
for the Stieltjes transforms, i.e. g(u) = l/(z — u), there is a constant c > depending only on V 
and g (one can choose c = 0(||g( 4 ) ||oo)) such that, for any N > 1, 



E„ 



g(Xi) - N p{u)g(u)du 



< c log N, 



(2.7) 



where fih is obtained by replacing V by V + in the definition of /i, and h is for example any 
iV-indepcndcnt smooth compactly supported function. We will now prove that this implies that 
(2.7) actually holds when replacing g by any smooth compactly supported 0, for example by a 
Helffcr-Sjostrand type argument, similar to Lemma 2.3. We can now apply formula (B.13) in [15] 



for the signed measure p 
correlation function of ph- We obtain 



-WW.) 

Pi 



p, with Stieltjes transform S, where p^ 1 ^^ is the one-point 



0(A)p(A)dA 



<C 



y<j>"{x)x{y)^S(x + iy)dx&y (2.8) 
/ (|^)| + |y||^(x)|) \ X '(y)\ \S(x + iy)\dxdy, (2.9) 



for some universal C > 0, and where x is a smooth cutoff function with support in [—1,1], with 
x(y) = 1 for |y| < 1/2 and with bounded derivatives. Note that x' is supported on 1/2 < \y\ < 1 
and cf), ((>' on compact sets, and that S is uniformly O ^ 1 °% N J on this compact integration domain, 

by (2.7), so the term (2.9) is easily bounded by O(||0|joo + ll^llcc) °ff N ■ Concerning the term 
(2.8), an easy calculation yields the bound j^(y ImS) = 0(l/y), so integrating from 1 to y we 

get \y lmS(x + iy)\ = 0(| logj/|) °% N , which is integrable, so (2.8) is 0(log N/N) as well, finally 
proving that (2.7) holds when replacing g by 4>. 
Following now Lemma 1 in [25], consider 

Z N (t) = E M ( exp ( T7 -^ 7 \J2 0(Ai) - N [ p(s)<j)( S )d S 



\osN 



, 2 

Then obviously log Z N (t) > 0, so 



log Z N (t) 



log Z N {t)~ log Z N {Q) < \t\-\ogZ N {t) 

dt 



log TV 



Mt^i/ log N 



N 
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so using (2.7) we get that Z^{t) < e c '*', from which Lemma 2.1 easily follows. rj 
2.3 Analysis of the loop equation 

This section analyzes the loop equation (2.10) in the following Lemma 2.2. Its proof is very similar 
to [6] except that, instead of the logarithmic Sobolev inequality which was valid only for convex 
V, we will use Lemma 2.1. Furthermore, since the support of the restricted measure pS N ' K ^ has 
changed, the integration contours in (2.16) are chosen slightly differently from those in [6]. 

In the form presented here, we follow closely the proof in [25]. We now introduce some notations 
needed in the proof. 

• ton is the Stieltjes transform of p[ N \s)ds, evaluated at some z with Im(z) > 0, and to its 
limit: 

m N (z)=E, (l-f^-L-) = f J_ p W(t)di, m(z)= [ —p(t)dt- 



s(z) = -2r(z)y/(A- z)(B - z), where the square root is defined such that 



f(z) = \J{A — z){B — z) ~ z as z — > oo; 

• bj\[(z) is defined by 

z - t 

finally, c N (z) = - N ^k N {z) + ± - lj m' N (z), where 

/ n 



k N (z) = var M ^ 



Z — A; / 
\fc=l / 

Here the var of a complex random variable denotes var(X) = E(X 2 ) — E(X) 2 , i.e. without 
absolute value unlike the usual variance. Note that |var(X)| < ~E(\X — E(X)| 2 ). 

The loop equation (see [16,17,25] for various proofs) is 

(tojv — m) 2 + s(mjv — wi) + bjy = cm- (2-10) 

In the regime where |mjy — m\ is small, we can neglect the quadratic term. The term is the 
same order as |mjv — m\ and is difficult to treat. As observed in [1,25], for analytic V (hence 
analytic bisr), this term vanishes when we perform a contour integration. So we have roughly the 
relation 

(rnjy - m) ~ ^ var M (j^ J^^j > ( 2 - n ) 

where we dropped the less important error involving m' N (z)/N due to the extra 1/N factor. With 
no convexity assumption on V, the difficulty will be to estimate the above variance to immediately 
obtain an estimate on ttin — to; this is the reason why we will introduce a convexified version of 
the measure /i in the next Section 3. To quantify more precisely (2.11) we will use the following 
result, already proved in [6] for convex V. 
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Lemma 2.2 Let S > 0. For z = E + irj with A + 5<E<B — 5 assume that 

-Lk N (z)^0 (2.12) 

as N — > oo uniformly in r\ > 7V _1+a /or some < a < 1. Then there are constants c, K > smc/i 
ttat /or any 7V- 1+a < -n < k, A + S < E < B - 5, 

\m N (z)-m(z)\<c(-^ + j^k N (z)y (2.13) 

Proof. First, for technical contour integration reasons, it will be easier to consider the measure 
(2.4) instead of here. More precisely, define 

kP(z) = var M (iv, K ) > 



Then it is a direct consequence of (2.5) and (2.6) that for any k > there is a constant c > such 
that uniformly on n > N~ 10 (or any power of N), 

|m^-m w | = 0(e- cAr ), \k%> ~ k N \ = 0(e~ cN ). (2.14) 

From now, we choose a fixed K > such that all the zeros of r are at distance at least 10k 
from [A, B] (this is possible because V is regular). Consider the rectangle with vertices B + 5k + 
LZV~ 10 , A - 5k + iiV~ 10 , A -5k- iiV~ 10 , B + 5k - iN~ w , call C the corresponding clockwise closed 
contour and £ the one consisting only in the horizontal pieces, with the same orientation. From 
(2.10), we obviously have, for z £ £', 

1 f (m N (Q - m(Q) 2 + s(Q(m N (Q - m( Q) + MO ~ cjv(Q ~ = Q 



2m J c , K0(«-0 

Note that the above expression makes sense for large enough TV, because then r has no zero on C. 
Using (2.14), this implies, for rj > N^ 1 , 



J_ f rfCj) - M£)) 2 + <0(mg(0 - m(Q) + MQ - 4%) _ n , 

2^ K0(*~0 4 ^ ( y 



Now, as p[ N,hlS> and p are supported on [A — k, B + k], — m and cjy are uniformly 0(1) in the 
vertical segments of C Consequently, from the above equation 

J_ f (mjg(|) - ^(0) 2 + 3(0(^(0 ~ ™(0) + MQ - jg(Q = 10 

As 6jv and r are analytic inside £, for z outside C we get 

1 f (m&>(0 - m(0) 2 + s(0rf(0 - m(Q)-4 K) (Q dc = o(JV _ 10) _ 



27ri7 £ r(0(«-0 
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Remember we define f(z) = y/(A — z)(B — z) uniquely by f(z) ~ z as z — > oo. Moreover, 



,(«) 



m|(z) = 0(2 2 ) as \z\ — > 00 because p and /9[ JV ' K ' are compactly supported 



\m}ff {z) - m(z)\ 



B + K 



A — K 



p(t) -/>(".»>(*) 



at 



B + K 



1 



(If 



(z- 2 ) 



Consequently, the function s(m^ — m)/r = — 2f(m { £ > — m) is 0(2; x ) as |z| — > 00. Moreover, it 
is analytic outside £, so the Cauchy integral formula yields 



2mj c 



r(0(*"0 



d{ = -2/(z)(< ) -m)(«) l 



proving 



2/(z)(m^(z)-m(z)) = -^ 



("*^(0-m(0) 2 -cW(0 



d£ + 0(7V~ 



K0(*-0 

Consider now the following rectangular contours, defined by their vertices: 

£1 : £> + 3k + i3k, A — 3k + i3K, A — 3k — i3K, B + 3k — i3K, 
£2 : B + 4k + i4K, A - 4k + i4«, A - 4k - i4«, 5 + 4k - i4«. 



(2.15) 



(2.16) 



In particular, note that all the zeros of r are strictly outside £2. For z inside £2 and lm(z) > N 1 , 
by the Cauchy formula, equation (2.15) implies that 

— 2s(z)(m£ (z) — m{z)) 
= -(m^(z)-m(z)f + c^ 



2?ri 



-dC + 0(iV- llJ ). (2.17) 



In the above expression, if now z is on £1, \z — £| > k, and on £2 |r| is separated away from zero 
by a positive universal constant. Moreover, eff (£) can be bounded in the following way. For any 
£ € £2, there is a smooth function supported on [A — 2k, B + 2k] which coincides with on 
[A — k, + k], Moreover, this choice can be made such that ||ff{||oo, Hfl^Hoo, Hs^'Hoo are uniformly 
bounded in £ € £2. Then 



1 



N 



var^(N,«) 



1 

iV2 



var^tiv,,) ^Xj.9?(Afe)^j 



1 



N 



\k=l 



(l + o(l)), 



where the last equality follows from (2.5). Now, from Lemma 2.1, this last variance is uniformly 
bounded by c (log N) 2 , with c uniformly bounded in £. This proves that k%\t) is 0((logiV) 2 /jV 2 ), 
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uniformly on the contour £2- Moreover, -4to)^ = 0(N 1 ), so finally c^(£) is uniformly 0(N 1 ) 
on £2 and (2.17) implies 

-2s(z)(m%\z) - m(z)) = -(m ( «\z) - m{z)) 2 {z) + O (sup - m| 2 ) + 0(iV" 1 ). 



Moreover, from the maximum principle for analytic functions, sup £2 
so the previous equation implies 



■m| < sup £i |m 



(«) 

AT 



sup |mj^ — m| = O (sup \rn$ — m\ 2 H 

£1 V £1 N 



We know that p[ N \s)ds converges weakly to p(s)ds (see [2]), so by (2.5) and (2.6) p[ N ' K \s)ds 
converges weakly to p(s)ds. On C±, z is at distance at least k from the support of both p[ N ' K \s)ds 
and p(s)ds so, on Ci, — m converges uniformly to 0. Together with the above equation, this 
implies that 

» 



sup \my - m\ = O 
£1 V v 

By the maximum principle the same estimate holds outside Ci, in particular on £2, so equation 
(2.17) implies that for z inside £1 



- 2a(z)(m%\z) - m(z)) = -(m%>(z) - m(z)f + c%>(z) + O 



(«)/ 



(2.18) 



Moreover, 



^l TO w ) ( z ) 
1 



1 

~N2 



3 



< — Im mft\z) < J^\™%\z) - m{z)\ + ^|Im m(z)| < ^|m^(z) - m(z)| + ^ (2.19) 

for some constant c. We used the well-known fact that Im m is uniformly bounded on the upper 
half plane 3 . On the set A + 8<E<B — 8 and \rj\ < «, we have inf \s\ > 0. Therefore (2.18) takes 
the form 



1 



1 



1 



(2.20) 



From the hypothesis (2.12), if iV- 1+a < 77 < k and A + <5<£;<B-J, then 



— m I < c\m 



(«) 

TV 



(2.21) 



for some c > and £jv 

,(«) 



as JV -> 00. For large N, (2.21) implies that |to (a 



N 



m\ < 2en or 

m\ > 1/c— 2ejv- Together with |to^? — m|(i£ + ire) — > and the continuity of |to^ — m\ in 
the upper half plane, this implies that \mf^ — m\ < 2em and therefore — m\ — > uniformly 



3 This follows for example from properties of the Cauchy operator, see p 183 in [9]. 
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on N 1+a < r/ < n, A + 5<E<B — 5. Consequently, using (2.20), this proves that there is a 
constant c > such that for any 77 > iV~ 1+a , A + 5 < E < B - 5, 



\m^(z)-m{z)\<c^— + 1 ^k^(z) 

The same conclusion remains when substituting (rcsp. kff) by (rcsp. kpj) thanks to 

(2.5) and (2.6). n 

To prove rigidity results for p,, the above Lemma 2.2 will be combined with the following 
Hclffcr-Sjostrand estimate, already proved in the following form in [6]. 

Lemma 2.3 Let 5 < (B—A)/2 and E e [A+5,B—S] andO < ?y < 5/2. Define a function f = f E . v : 
R — > M such that f(x) = 1 for x £ (— 00, E — rj\, f(x) vanishes for x G [E + rj,oo), moreover 
\f'( x )\ < o] 1 an d \ f"{ x )\ ^ crj~ 2 , for some constant c. Let p be an arbitrary signed measure and 
let S(z) = f(z — x)~ 1 p(x)dx be its Stieltjes transform. Assume that, for any x £ [A + 5/2, B — 5/2], 

\S(x + iy)\ < — for 77 < y < 1, and | Im S(x + iy)\ < — for < y < r). (2.22) 
JSy l\y 

Assume moreover that J R p(X)d\ = and that there is a real constant T such that 

|Ap(A)|dA < (2.23) 

Then for some constant C > 0, independent of N and E € [A + 5, B — 5] , we have 

CU\\og V \ 



f E (X)p(X)dX 



< 



N 



3 Convexification 

3.1 Outline of the main ideas 

The Hamiltonian H = Hn of the measure \x ~ exp(— (3 NO.) is given by 

1 N 1 

fe=l \<i<j<N 

TL is not convex, but its Hessian is bounded from below, V 2- H > — W. We will modify this 
Hamiltonian by an additional term 

e n 
H:=H + (W + 1)^X1 X a = N- 1 / 2 J2(9 a (X j )-g a (y j )), (3.1) 

Q=l j=l 

where the real functions g ai a = 1, 2, . . . ,£, will be determined later and will be independent of N. 
Here we denoted by a slightly modified version of the classical location of the points, defined by 
the relation 

£ p(s)ds = 3 -j^, j = l,2,...,N. (3.2) 
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Compared with 7^ defined in (1.6), there is a small shift in the definition which makes a technical 
step (Lemma 3.3) easier in this section. In all estimates involving jj this small shift plays no role 
since maxj |7j < CN~ 2 / 3 . In particular the crude large deviation bound (1.7) holds for 7's 

as well: 

P„ (3* 6 [1, N} I |A fe - 7*| > e) < Cl e^ N . (3.3) 
The TV -1 / 2 normalization in the definition of X a is chosen such that the vector 
G Q :=^V- 1/2 ( 5 o(7i),ffo(72),...,ffo(7iv)) tR N 

is £ 2 -normalized. 

Define the random variables 

A := max U £ | A, - 7, |, ± J> " ^ ) * AW : = 5 + 7^ " ^) 2 



for any < <5 < 1. Clearly 



|^|<iV 1/2 ||.9LllooA. 



We then have, for any vector v € R , that 



<v,(V 2 H)v) = l£ 



1 

2(iy + l)]T 



> 



-e 



(Ui-Vj) 



iV ^ (A* - Xj) 2 
i<j ■>> 



(W + l)J2\(G a ,v)\ 2 



^ + 2A(M/ + l)E(ll.9«l| 2 oo + ll.^l|oc:||.^||oo) 



(3.4) 



where we used a simple Schwarz inequality 



i j 3 

> (e^k) 2 - 2 E [^(Ai) -g' a m 2 \\v\ 



in the last step. 

We will define below a nonncgative symmetric operator Q on C N via a quadratic form 

jv 

(v, Qv) = E %(«i-^) 2 



(3.5) 
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such that for typical point configuration A = (Ai, A2, . . . , Xn) we have 



N 



-■J 1 



> 



^ Qij{vi - vj) 2 . (3.6) 



N (A, - A,-) 2 

In our applications, we will then choose £ to be a large but JV- independent number, we will let G a , 
a = 1, 2, . . . , I be the eigenfunctions corresponding to the lowest I eigenvalues Hi < /12 < ■ ■ ■ < He 
of the nonnegative operator Q. Thus we will have the operator inequality 

i 

Q + (W + 1) \G a )(G a \ > mm{(W + 1), 

a=l 

In Section 3.2 we will show that for I sufficiently large, independent of N, we have ne+i > W + 1. 
Setting 

1 

C(l) := 2(W + 1) ^ (II^HL + Rllool^lloc), 

a=l 

we will obtain from (3.4) that 

(v, (V 2 H)v) >({W + l)-W- C(*)AW) ||v|| 2 > (3.7) 

as long as < 1/C(£). From now, we choose 5 = 1/(2C(^)) and as TV — > 00, we have < 
1/C(£) with very high probability, thanks to the large deviation estimates (3.3). 

To prove that ^e+i > W + 1, we only need to estimate the low lying eigenvalues of Q and 
we need to understand the low lying eigenfunctions G Q . Since the only requirement for Qij is to 
satisfy the bound (3.6), we have a substantial freedom in choosing Qi j conveniently. There are 
many ways to choose Qij ; we will give one possible approach that relics on enlarging the space by 
a reflection principle in the next section. Roughly speaking, we will construct an operator 1Z with 
periodic boundary conditions on the set consisting of the original set and its "reflection" . Wc then 
choose Q to be the restriction of 1Z to the symmetric (under the reflection) sector. The operator 1Z 
is translation invariant, hence it can be diagonalized via Fourier transform and the eigenfunctions 
are explicit. The reader may skip the next section on first reading as it contains fairly elementary 
arguments that are independent of the rest of the paper. 

3.2 Slow modes analysis 

Let / := [1, NJ be the index set of the vectors v. The original operator Q is defined on the space 
£ 2 (I). We enlarge this space to £ 2 (I), where / := [— N + 1,NJ. We extend any vector v € £ 2 {I) 
by reflection to a vector v e £ 2 {I) as follows 

Vj=Vj, j = l,2,...,N; (3.8) 
Vj=V!-j, j = 0, ...,-N+l. (3.9) 

We will often view the set / modulo 27V periodic, i.e. we consider it as 2iV points on a circle and 
identify —N with N. We can thus also view £ 2 {I) as the space of vectors with periodic boundary 
condition V-m = vn- The algebraic operations on the indices will be considered modulo 27V. 

We consider the natural translation invariant distance on /. Define the function m(n) for neZ 
such that m(n) E ^—N + 1,NJ and m(n) = n mod(27V). Then the distance between k,£ e / is 
defined as d(k,£) := \m(k — £)\ which ranges from to N. 
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Lemma 3.1 Let e > be sufficiently small, depending only on V. Define 

R ki= R kA~jj d{k t! 3 2 » k,eeI=[-N + i,N}, (3.10) 

and 

Qi e J ( l ■ ■■= ^ + *>' u + + i [i, Nl (3.ii) 

Then there is a constant c\ > depending only on V such that for any e > there is a constant 
C2 > (depending on V and e) such that for any N and i,j G [1, NJ 

The relation between Qij and Rij is dictated by the requirement that 

(v, Qv) P{I) = {v,Kv) p{J) := R iAvi - «i) 2 : Vv e £ 2 (I) 7 

which can be easily checked from (3.8). 

Proof. Recall that [A, B] is the support of p, p > on (A, B) and p has a square-root singularity 
at the two endpoints, i.e. it vanishes as p{x) ~ saV x — A as x — > A + and p(x) ~ sb\/B — x as 
x — s- B~ with some positive s Ai sb- 

From the large deviations of the extreme eigenvalues (included in (1.7)), we know that for any 
K > there is a c(k) > such that 

Ai < A - k) + P p (\ N > B + k) < e- c(K)N . (3.12) 



Fix a positive number s < min(s^i, Sb). Then there is a kq > 0, depending only on V, such 
that 



p(x) > s\Jx - A - t xe[A , A+tio] + l xe [ A+KO<B _ Ko ] +sVB-x-1 X £[b- KOi a]- (3-13) 

Let e < ck,q^ 2 with a small positive constant c. Suppose that k < iV/2; if fc is near the upper edge, 
the argument is similar. Since 

™ fc - 5 /" A * 1 nt^ , 1 



2 

P = 



N^2 6x ™ ON' 



from the first relation we get 

c(fc/iV) 2/3 < % - A < C{k/N) 2 ^ 3 (3.14) 

with some positive constants c, C, depending only on V. Subtracting the first and second relations 
and using (3.3), we obtain that for any fixed K 



P 



< x (3.15) 
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apart from an event of exponentially small probability (i.e. of type exp(— c(e / K)N)). 

Additionally, assume now that k > Ne. Under (3.15) we easily see that Xk £ (7fc/2j 73fc/2)> 
since both f~* 2 P and J^ k/2 P arc of the order k/N which is larger than e/K if K is large enough 
(depending only on V). Then (3.13) and (3.15) imply 

|A fc - 7fc| < Ce(% - A)- 1 ' 2 < Ce(k/N)-^ 3 , Ne < k < N/2 

with exponentially high probability and with a constant C depending only on V. 

Now we consider the k < Ne case. Using (3.12) with k = e 2 ! 3 and (3.14), we have (apart from 
an event of exponentially small probability) 

% - Xk < C{k/N) 2 ' 3 + k< Ce 2 ' 3 . 

Finally, still when k < Ne, i.e. % < A + Ce 2/3 then (3.15) implies that X k < A + de 2/3 with a 
large C\, i.e. 

X k -%<Ce 2 / 3 , k<Ne, 
still apart from an event of exponentially small probability. Summarizing all cases, we obtain that 

holds with overwhelming probability. 
Now let |i — j'l > Ne, then 



l*-Jl 
AT 



1i 



P 



>c|7i-7j|e 1/3 



since either i or j is larger than Ne and smaller than N(l — e), say Ne < i < N(l — e), and then p 
is at least of order e 1 / 3 in the neighborhood of 7V If \i — j\ < Ne, then we have the trivial bound 
\li < Ce 2 ! 3 . Combining these, 

holds for any Furthermore, clearly \Xi — 7z| < Ce 2 ! 3 from (3.16), so we have proved that 

with overwhelming probability and for any i,j. In other words, there is a constant C (depending 
only on V) such that for any sufficiently small e and for some c(e) > we have for any N and 

ye n,m 

( 1 rv 2 / 3 \ 

< < e- c(e)N - (3-17) 



The proof of Lemma 3.1 will therefore be complete if we can prove that 

\i- i\ < min{d(i,j),d(l - i,j),d{i,l-j),d(l - j)}, i,j £ I. 
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Note that we can assume i > j and then \i — j\ = \m(i — j)\ = d(i,j) is obvious. Moreover, as d is 
symmetric to the reflection i — > 1 — i on /, i.e. d(k, i) = d(l — k, 1 — £), we just need to prove that 
< d(l - i,j). If i + j < N + 1, then d(l - i, j) = i + j - 1 > \i - j\. If N + 2 < i + j < 2N, 
then d(l - i, j) = \m(i + j - 1)| = \i + j - 1 - 2N\ =2N+l-i-j>i-j because i < N, 
completing the proof. rj 

We use the matrix Q in the previous lemma instead of bounds of type (3.17) because it is 
related to R, a circulant matrix, allowing to derive its eigenvalues and eigenvectors in an explicit 
way. 

Lemma 3.2 Let R = R( s > be the matrix (Ri,j)u where the matrix elements are defined 

in (3.10). Then the eigenvalues v\, . . . , V2N of R are 

2N-1 

and the associated normalized eigenvector is w£ = (2-/V) -1 / 2 (e l27r:, 2^)j = _jv+i,...,Ar. 

In particular, for any given W > there is a sufficiently small e such that for large enough N 
we have V2N > W. Moreover, for any given e > and s > there is some a > depending only 
on e and s such that for any N 

{k : |i/ fc | > s} C [1, oj U [2JV + 1 — a, 2Nj. 

We remark that the matrix R defines a symmetric operator w — > i?w and a quadratic form 
(w, Rw) = Ri,jWiWj in £ 2 (I). It is related to the quadratic form TZ via 

(w,ftw) = ^R itj {wl + w 2 ) 2^ /.•, ,"•,"•., = 2^ 2A r||w|| 2 - 2(w,i?w), (3.18) 

i.e. we have 1Z = 2&2N — 2i? on ( 2 (I). 

Proof. The first assertions, about the eigenvalues and eigenvectors, is a general fact about 
circulant matrices and can be obtained by Fourier transform in |0, jjy, . . . , 2 ^ 1 }. 
Concerning the distribution of eigenvalues, note that 

1 e 2 / 3 2/3 f 1 dx _ 1/3 f 1 '* _dx_ 

We therefore have, for sufficiently small e > 0, V2N > W for large enough N. We now write 

r 2/3 



e=o n 2 



0<j<2N- 1, 



with the convention 6 2 Af = 0. Then note that | a, | < 2|1 — e l27r 2N | 1 f and write ^ = Xw=o 1 a j(bj~ 
bj+i) to obtain 



JV|1 -e i27r w| 



e 2 + (^ N2 



iv/ V ./V 



< c(e) 



iV|l - 
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Consequently, if \i/ k \ > s then |1 — e l2lT 2» | < c(e,s)/N, which only happens for indices k whose 
distance to or 2N is uniformly bounded. q 

Lemma 3.3 For a fixed at N, there exists an N -independent function g a with \\g a \\oo + ||<?LII°° + 
||.9q||oo < oo such that 



Define 



g' a (%) = V2co S ^(A-^Q. 
G a = VV- 1/2 ^(7i), • ■ • = V27V- 1 / 2 (cos (n(2j - 1)^)) . =j 



it is easy to see that ||G Q || = 1. For any M > there is e > and £ > such that, for large 
enough N we have, as an inequality between positive operators on £ 2 {I), 

i 

Q + AlJ2\ G a )(G a \>M, 

a=l 

where Q = was defined in (3.5) with coefficients Qij = defined in (3.11). 

Proof. The existence of g a follows easily from the fact that the density p has a square- 
root singularity near the edges, see (1.4). In fact, we can choose the functions g a such that 
9a( x ) = a/2 cos (^irsACt(x — A)^ 2 ) for x near A, and a similar relation holds near the upper edge. 
As (v,ftv) = (v, Qv), ||vj| 2 = 2j|v|| 2 and 

|(v,G tt )| 2 = |(v, UQ )| 2 = Kv, M2W _ Q )| 2 , 

we just need to prove that the operator inequality 

2K + M \u k )(u k \>M 

fee[l,^]Ul2JV-^,2JV-l] 

holds in £ 2 (I) for some fixed constant £ and for large enough N. This is equivalent to the statement 

M M x , 

R + — <V2N + — 22 l«*)(«fc|. (3.19) 

fce[l/]U[2iV-^,2iV"-l] 

since 1Z = 2v2N — 2R, see (3.18). Recalling that v k are the eigenvalues of R — R^ e \ we need to 
check that for some fixed e > and £ we have 

M M 

v k + — < V2N + -|-lfee[l/]U[2iV-^,2iV-l] 

for any k. Since V2N is the top eigenvalue, this inequality is obvious if k £ [1,^] U f2N — £, 2N — 1]. 
Moreover, Lemma 3.2 proves the existence of some fixed e > and £, such that v^n > M and 
{k :\v k \< |M} C \l, Ij U {2N - t,2N - 1] hold for large enough N (depending on M as well as 
e and £). This concludes the proof. rj 
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3.3 The locally constrained measures 

In this section some arbitrary e, a > are fixed. Let 8 be a continuous nonncgativc function with 
6 = Oon [-1, 1] and6>" > 1 for \x\ > 1. We can take for example 6{x) = (x-l) 2 l x> i + (a;+l) 2 l a; <-i 
in the following. 

Definition 3.4 For any s,£> 0, define the probability measure 

Av {s,l) = e -PNH v ._ _}_ e -fSN^-pNJ2 t y^-pN(W + l)j:Uxl d ^ (3 2Q) 

where 

• the function g a was defined in Lemma 3.3; 

. ^(A) = ^(^E l =i(A,-7 i ) 2 ); 

• ipi t j(X) = jj9 (\fc~x N Qi y j(Xi — Aj)), where c\ was defined in Lemma 3.1. 

Lemma 3.5 There are constants c,£,s > depending only on V, such that for N large enough 
i/( s '^) satisfies, for any v € R N , 

(v,(V 2 H„)v) >c||v|| 2 . 

Proof. Using the notation (3.1), we have % v = H + iftW j VVj U P to an additive constant, 

so 

(v,(V 2 H„)v) =(v(V 2 7i),v) +Cl J2Qi,jO" (\ r, A (>,.,■ X, - A,)) (v z - v 3 ) 2 

+^'(£X>-*> 2 ) (e^-^) +2^'^E(A^70 2 )iivii 2 . 

(3.21) 

We now use (3.4) to get a lower bound for (v(V 2- H),v). Note that the second 9" term is positive 
and that for x > 0, 6'(x) > 2(x — l)l x >i; we therefore get the following lower estimate of (3.21): 

E ( (A . r A )2 + c i ^ Q'.j 1 ( (A . 1 A )2 < c i ^ Qm)) fa -"j) 2 +( ff + i )Ei( G ^ v )r 

2, J Of — 1 

- (V + <7(£)A«) ||v|| 2 + 4s (sA - l)l sA _ 1>0 j|v|| 2 , 

which is greater than 

£ 

(v, (ciQ + (W + 1) XI |G a ><G a | - (W + l))v) + (1 - C(£)A^)\\v\\ 2 + As(sA - l)l sA _ 1>0 ||v|| 2 . 

a=l 

Choosing M — (W + l)/ci in Lemma 3.3, for e small enough and £ large enough the above scalar 
product term is positive for any v and large enough N (note that c\ does not depend on e) . For this 
choice of £, taking s = C(£) makes the other terms all together positive (without loss of generality 
we can assume C(£) > 1/4), concluding the proof. rj 
From now, we abbreviate v for z/'w, where s and £ are fixed such that the conclusion of Lemma 
3.5 holds. 



18 



3.4 Equivalence of the measures v and /i 

Wc say that a sequence of events (An)n>\ is exponentially small for a sequence of probability 
measures (mjv)jv>i if there are constants 5,ci,C2 > such that for any TV we have 

m N (A N ) < Cl e- C2 < 

Lemma 3.6 For any s,£ > 0, the measures (/-i^ n,> )n>i an d (^ s '^' N ')n>i have the same exponen- 
tially small events. 

Proof. First note that T-L v > H^, so Z v < Z^. We claim that the following inequality holds: 

hgZ li <logZ v + Er(pN(Hv-'H li )). 
To prove it, by Jensen's inequality we have 

l0g J 6 "^ Je^dX ^- J WH'-KJ je-WH^ - 
We now bound 

K M (N(H V - H M )) = M NX l) + ^(%) + M N ^ (S) ) ( 3 - 22 ) 

in the following way. By Lemma 2.1, E^AA 2 ) < Co,(log A) 2 ; by Lemma 3.1 (together to (2.6)), 
Efj,(Nipij) is subcxponentially small for any indices i and j; finally E M (A?// S )) is also subexponen- 
tially small by (3.3) (together to (2.6)). Altogether, we get that there is a constant c > such that 
for any N > 2 

log Z v < log < log Z v + c(log A) 2 . (3.23) 
Let (Ajv)at>i be now a sequence of events exponentially small for /x. By (3.23) we have 

K(A N ) < e« l °z N )\(A N ), 

so (An)n>i is also exponentially small for v. 

Assume now that (Atv)at>i is exponentially small for v. there are constants 8,c\,C2 > such 
that for any TV we have 

Then for any t we have 

P^(A N ) = F^A N n {0N(H v - > t}) + Pp(A N n {0N(H„ - H M ) < t}) 
< ¥^{{PN{U V - Uy) > t}) + e^ v (A N ), 

where we used Z v < Z^. Choosing t = A" 5 / 2 makes the second term exponentially small, and the 
first one as well by using as previously Lemma 2.1, Lemma 3.1 and (3.3). rj 

From the previously proved equivalence of the measures [i and v 1 we can easily obtain rigidity 
of the particles at scale A -1 / 2 , 

Proposition 3.7 For any a,e > 0, there are constants S,ci,C2 > such that for any N > 1 and 
k g laN,(l-a)Nj, 

P»(\\ k -Jk\>N- 1 i+^<c 1 e- c * NS . 
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Proof. From Lemma 3.5 about the convexity of T-L U: we get by the classical Bakry- Emery 
criterion [3] that v satisfies a logarithmic Sobolev inequality with constant of order 1/2V, so by 
Herbst's lemma concentration at scale TV -1 / 2 holds for individual particles for v. there is a constant 
c> such that for any N > 1, k € [l,iV] and x > 0, 

P,(|A fc -E, y (A fc )| >xN~^ <e~ cx . 

By Lemma 3.6, this implies that for some constants <5, ci,C2 > 0, 

P M (|A fc -E„(A fc )| >A^+ £ ) <c ie - C2A,i . (3.24) 

From the above equation we get |E„(Afc) — E jU (Afe)| is of order at most A r_1 / 2+£ , so (3.24) holds 
when replacing E„(Afc) by E M (Afc), proving concentration at scale TV -1 / 2 for \i. 



Define 7^ by 



(JV) 

^ P{ N) = ± (3-25) 



The proof will be complete if we can prove that for any e > and k £ [aiV, (1 — a)7Vj, for large 
enough N we have 

| 7 f ) - 7fe |<^ 1 / 2+e , (3.26) 

By Lemma 2.2, |mjy — m| — > for 77 > iV _1 / 2+£ , because on this domain -jbkpf — > 0, as concen- 
tration at scale JV -1 / 2 holds for fi. So using Lemma 2.3 we finally get that (3.26) holds, finishing 
the proof. rj 



4 The multiscale analysis 

The purpose of this paragraph is to prove the following proposition: if rigidity holds at scale 
N~ 1+a , it holds also at scale N~ 1+ 4 a . The argument very closely follows Section 3.3 of [6] and we 
will just explain the modifications. 

Proposition 4.1 Assume that for some a € (0, 1) the following property holds: for any a,e > 0, 
there are constants 5, c\, C2 > such that for any N > 1 and k € [aiV, (1 — a)-/V], 

(|A fc - 7fc| > A^ 1+a+£ ) < c ie -^- NS . (4.1) 

Then the same property holds also replacing a by 3a/4; for any a,s > 0, there are constants 
<5, ci, C2 > such that for any N > 1 and k € [ccJV, (1 — a)iV], we have 

P M (|A* - 7*1 > N- 1+ i a+ ^) < Cl e~^ N \ 

Proof of Theorem 1.1. This is an immediate consequence of the initial estimate, Proposition 
3.7, and iterations of Proposition 4.1. rj 

As in Section 3.3 of [6], two steps are required in the proof of the above Proposition 4.1. First 
we will prove that concentration holds at the smaller scale iV~ 1+ t. 
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Proposition 4.2 Assume that (4-1) holds. Then for any a > and e > 0, there are constants 
ci, C2, <5 > sitc/i £/ia£ /or any TV > 1 and fc €E [aiV, (1 — a)./V], 



|A fc -E M (A,)|>^) 



After the better concentration from this proposition, the rigidity can be improved to the scale 

Proposition 4.3 Assume that (4-1) holds. Then for any a > and e > 0, there is a constant 
c > such that for any N > 1 and fc G [aJV, (1 — a)iV], 



(AT) 

7 fc -7fc 



< c- 



JV 



where jj. is defined in (3.25). 



Propositions 4.2 and 4.3 are the equivalent versions of Propositions 3.12 and 3.13 of [6] with no 
convexity assumption on V. Proposition 4.1 can be proved exactly in the same way as Proposition 
3.11 [6] by using Propositions 4.2 and 4.3. Notice that this argument does not use the convexity 
of V. We now explain the proof of Propositions 4.2 and 4.3. 

The convexity of V is used critically in the proof of Proposition 3.12 of [6]. Our measure \i is 
not convex, but thanks to Lemma 3.6, it has the same exponentially small events as the measure 
v {s,£) £ or an y fi xec j S ^g > o. Hence it suffices to prove (4.2) with [i replaced by v^ s ' e \ Choose an 
appropriate s, £ such that the Hamiltonian of v = is convex (Lemma 3.5). Then the proof of 

(4.2) with fi replaced by v is identical to the proof of Proposition 3.12 of [6] since the measure v 
is now convex. 

For the proof of Proposition 4.3, we can follow the proof of Proposition 3.13 in [6] line by line. 
At a single place, in estimating the second term on the r.h.s. of (3.51) in [6], the spectral gap 
inequality for (Eq. (3.12) in [6]) was used, but the necessary estimate immediately follows from 
Proposition 3.7. 
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